Quantum phase transitions and quantum fidelity 
in free fermion graphs 
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In this paper we analyze the ground state phase diagram of a class of fermionic Hamiltonians by 
looking at the fidelity of ground states corresponding to slightly different Hamiltonian parameters. 
The Hamiltonians under investigation can be considered as the variable range generalization of the 
fermionic Hamiltonian obtained by the Jordan- Wigner transformation of the XY spin-chain in a 
transverse magnetic field. Under periodic boundary conditions, the matrices of the problem become 
circulant and the models are exactly solvable. Their free-ends counterparts are instead analyzed 
numerically. In particular, we focus on the long range model corresponding to a fully connected 
directed graph, providing asymptotic results in the thermodynamic limit, as well as the finite- 
size scaling analysis of the second order quantum phase transitions of the system. A strict relation 
between fidelity and single particle spectrum is demonstrated, and a peculiar gapful transition due to 
the long range nature of the coupling is found. A comparison between fidelity and another transition 
marker borrowed from quantum information i.e., single site entanglement, is also considered. 



I. INTRODUCTION 

In the last years, several papers have shown that the 
study of quantum phase transitions (QPTs) can benefit 
from the application of concepts borrowed from quan- 
tum information theory. In particular, the behavior of 
entanglement at quantum critical points, due to its con- 
nection with quantum correlations, typically presents pe- 
culiar features which can be used as transition markers, 
thereby providing further insight into these phenomena 
P|. More recently, QPTs have also been analyzed from 
the point of view of Berry's geometric phase |2(, a quan- 
tity which has in fact been found to contribute to entan- 
glement entropy 

One of the most recent reasons of interest in quan- 
tum phase transitions for the quantum information com- 
munity is the relationship between criticality and deco- 
herence Another relevant aspect is connected with 
adiabatic quantum computing [5j, as the proximity to 
QPTs could cause the failure of adiabaticity conditions. 
On the other hand, crossing of QPTs could be desirable 
to transform easy-to-prepare states into non-trivial states 
corresponding to the solution of computationally relevant 
problems. Transitions with non-vanishing energy gap, 
like topological QPTs, could allow such an implementa- 
tion of adiabatic algorithms and therefore the achieve- 
ment of significant computational speedups 0. 

In Refs. 0, II] an approach to quantum phase transi- 
tions in terms of ground state fidelity (the overlap mod- 
ulus) was proposed. In practice, this consists in analyz- 
ing the fidelity between ground states corresponding to 
slightly different Hamiltonian parameters. The intuitive 
idea underlying this method is that in the proximity of 
critical regions a small change of the Hamiltonian param- 
eters can give rise to dramatic ground state (GS) varia- 
tions, due to the strong difference of the GS structure in 
opposite phases. The overlap between neighboring states 



is then expected to decrease abruptly at the phase bound- 
aries. In Ref. this idea was successfully applied to the 
study of the Dicke and XY models, where a one to one 
correspondence was found between the sudden drops of 
the overlap function and the critical points of the phase 
diagram. In Ref. || a more extensive analysis on general 
free Fermi systems was presented, providing a descrip- 
tion of the fidelity features in terms of the matrices of 
coupling constant entering the system Hamiltonian. In 
particular, explicit connections between the vanishing of 
the lowest single particle energy (the gap) and the fidelity 
drops were shown. An exemplification of these concepts 
in a model corresponding to a totally connected fermionic 
graph was also included. 

In this article we further develop the general formal- 
ism presented in Ref. , deriving in detail and expanding 
some of the results reported there. In particular, we char- 
acterize the fidelity behavior by a proper combination h 
of its second derivatives, namely the lowest eigenvalue of 
the Hessian matrix. Furthermore, we describe a class of 
exactly solvable models which offer the ideal playground 
for the study of the relations between ground state fi- 
delity and QPTs. The corresponding Hamiltonians are 
the cyclic variable range generalization of the fermionic 
Hamiltonian obtained by the Jordan- Wigner transforma- 
tion of the XY spin-chain in a transverse magnetic field. 
We also consider their free-ends counterparts, which are 
however analyzed numerically, as the analytical results 
available for the cyclic case make crucially use of the 
simple properties of the circulant matrices appearing in 
the presence of periodic boundary conditions. We then 
focus on the long range model corresponding to a fully 
connected fermionic graph. Along with some first order 
QPTs, two second order QPTs are found, as anticipated 
in Ref. Q for the free-ends case. One of them is an 
example of a gapful second order QPT, an uncommon 
situation also found in topological QPTs. The existence 
of a transition in the absence of a vanishing single parti- 
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cle energy turns out here to be due to the (unphysical) 
long range nature of the coupling. In order to ascertain 
the second order nature of such transition, in the ana- 
lytical cyclic case we also study the ground state energy 
derivatives, finding explicitly a discontinuity at the phase 
boundary in the thermodynamic limit (TDL). Further- 
more, the finite-size scaling behavior typical of second 
order transitions is recovered in the function h. Finally, 
we calculate the single-site entanglement and compare its 
behavior with that of fidelity in critical regions. 

The paper is organized as follows. In Sec. II we re- 
view and discuss in more detail the general results on 
quadratic fermionic Hamiltonians already presented in 
Ref. • In particular, after a brief summary of the diag- 
onalization procedure of Ref. @, Subsec. Til Al is devoted 
to the general expression of the GS in these systems, 
Subsec. Hi Bl to the explicit calculation of the fidelity, and 
Subsec. Til Cl to fidelity derivatives and the function h. In 
Sec. Ill we introduce the class of variable range models 
we have developed. We begin by recalling the fermionic 
Hamiltonian obtained by the Jordan- Wigner transfor- 
mation of the XY spin-chain in a transverse magnetic 
field, which motivated the introduction of our models, 
and then provide the analytical solution for the variable 
range cyclic case. In Sec. IV the long range fully coor- 
dinated system is described in detail, first providing an 
exhaustive analytic description of the cyclic graph (Sub- 
sec. IIV A|) and then presenting the numerical results ob- 
tained for the free-ends one ( Subsec. IIV B|l . In this sec- 
tion all the formalism developed in Sec. I for the fidelity 
approach is concretely applied to the study of the phase 
diagram of the system. The finite-size scaling of h in the 
neighborhood of second order QPTs and its asymptotic 
behavior in the thermodynamic limit are discussed, along 
with the GS energy. Then, in Sec. V we report our results 
on single-site entanglement, whose behavior in critical re- 
gions also shows transition signatures. Finally, we draw 
our conclusions in Sec. VI. Some technical details can be 
found in the Appendices. 

II. GROUND STATE FIDELITY IN 
QUADRATIC FERMIONIC HAMILTONIANS 

The most general quadratic Hamiltonian describing a 
system of L free fermionic modes is given by 

L L 

# = E 4^ + 2 E ( c t B v4 + h - c -) > w 

where Cj, cj are the annihilation and creation operators 
of the i-th mode and the matrices A 7 B are hermitian 
and skew-symmetric respectively, due to the hermiticity 
of H and the anticommutation properties of the fermionic 
operators. In the following, we restrict ourselves to the 
case where A and B are real matrices. The generality of 
the system described by Eq. Q can then be appreciated 
by representing A and B in terms of the generic (real) 



matrix Z, i.e., A = (Z T + Z)/2 and B = (Z T — Z)/2. 
This allows to identify the space of coupling constants 
of Eq. (JJJ with the L 2 -dimensional full matrix algebra 
Ml(JR). Correspondingly, as it will be illustrated below, 
the system's ground state properties and their changes in 
the presence of QPTs will be completely characterized by 
the matrix Z, which is hence of fundamental importance 
in our description. 

The above Hamiltonian can be rewritten in the more 
compact form H = (*tC* + TrA)/2, where * = 
(c 1; . . . , c L , c\, . . . , c\) T , = (c-l, . . . , c L , c 1; . . . , c L ), 

and C = a z <S> A + ia y ® B. This notation, also used 
in Ref. ficj , will be adopted throughout this section. 

The diagonalization problem for Eq. (JTJ has been 
solved in Ref. 0. This well known procedure is very 
powerful, as it involves operations in the parameter space 
of dimension L 2 , instead of in the full Hilbert space of 
the system, which has dimension 2 L . Below, we derive 
the results of Ref. by making use of the general prop- 
erties of real canonical transformations, which we briefly 
recall also in view of some simple applications in the next 
Subsection. 

Real canonical transformations - We consider the linear 
transformation 

L 

c 'i = J2( u jk c k + v jk c l) > ( 2 ) 
fc=i 

where u, v are real matrices. Correspondingly, '4'' = V\I r , 
where V — I2 <8> u + a x (£> v and the subscript is used 
to specify the dimension of the unit matrix whenever 
different from L. The transformation is canonical if it 
preserves the anticommutation relations. This implies V 
to be orthogonal, i.e., the canonical conditions 

uu T + vv T = 1 , (3a) 
uv T + {uv T ) T = (D . (3b) 

The Hamiltonian becomes H = (^C"*' + TtA)/2, 
where C = VCV T = a z ®A' + ia y <g> B'. Apart from the 
constant energy shift, which is still given by TxA/2, the 
structure of H is then preserved in terms of the primed 
quantities. One can thus define Z' — A' — B', therefore 
finding by straightforward calculations the simple rela- 
tion 

Z' = (u + v)Z{u-v) T . (4) 

Hamiltonian diagonalization - For the diagonalization of 
Eq. Q one then needs to diagonalize the 2L x 2L sym- 
metric matrix C, i.e., C = V T C'V with C" diagonal. We 
note that the eigenvalue matrix C" can always be writ- 
ten in the form C — a z ® A, where the L x L diagonal 
matrix A = diag(Ai, . . . , Ar) is positive semi-definite. 
Indeed, it is easy to see |2(j that det(C — ct.21) = 
imposes det(C + c^l) = as well, i.e., the eigenvalues 
of the symmetric matrix C appear in pairs of real num- 
bers of opposite sign. On the other hand, one also has 
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C = cr z <E> A' + io-y g> B', so that A' = A, B' = and 
hence Z' = A. From Eq. Q, by defining <I> = it + u and 
^ = u — v, one finally gets the important equation 



$Z* T = A , 



(5) 



also implying ^>Z T ^ T — A. Note that, due to the canon- 
ical conditions 10) obeyed by u and v, the matrices i> 
and "J must be orthogonal. It is then straightforward 
to derive the further relations A* = A$ = tyZ T 
and $ZZ T = A 2 <E>, ^Z T Z = A 2 *. In practice, one can 
for example solve $ZZ T = A 2 $ for <I> and A, and then 
obtain = A _1 <I>Z. Whenever Aj — 0, the equation 
Avp = <&Z cannot be completely solved for but one 
can recover the j-th. row of $ from A$ = ^Z T . Alter- 
natively, one first solves for \& and A, and then for $. 
Clearly, due to Eq. |J5J, <f> and \? are not independent, so 
that in general they cannot be calculated by diagonaliz- 
ing ZZ T and Z T Z separately. 

In conclusion, after solving Eq. © for <&, if? orthogonal 
and A diagonal, one introduces g = ($ + ^)/2 and h = 
(<f> — ^)/2, in terms of which the canonically transformed 
operators diagonalizing the Hamiltonian are defined by 
Vj = Z)fc=i(3ifc c fe + ^7fe c l)- The Hamiltonian finally reads 



h = J2a 



En 



(6) 



where Eq = Tr(A — A)/2 is the ground state energy and 
the Aj's are found to play the role of single particle en- 



ergies. 



A. Ground State 

In this subsection we present the explicit expression of 
the ground state of Eq. JIJ, starting from the Ansatz of 
Ref. The Hamiltonian, although not number con- 

serving due to the terms containing the antisymmetric 
matrix B, preserves the number parity Pm = (— 1) , 
where N = 53j=i c ] c j 1S the fermion number operator. 
In formulas, [H, N] ^ and [H,P N ] = 0. Hence, the 
ground state must be a superposition of states with even 
or odd number of fcrmions exclusively. In the case of 
even parity, it can be of the following form 



I*; 



- Afexp 




(J) 



where Af is the normalization factor, |0) is the fermionic 
vacuum (cj|0) = 0) and G is a real L x L antisymmetric 
matrix. The link with the above diagonalization proce- 
dure can be made by imposing the ground state condition 
r) k \^>z) = 0, V k H3. Then, the matrix G must be the 
solution of the equation 

gG + h = , (8) 



where it is worth recalling g = ($ + ^)/2 and h — 
(<t> — ^)/2. Note that the latter equation is not al- 
ways solvable, corresponding to cases where either the 
Ansatz {7J does not hold or the ground state parity is 
odd. We will return on this point below. 

The bridge with the picture of the previous subsection 
can be made more stringent whenever A (and hence Z) 
is invertible. In fact, by using Eq. (J5J one can write g = 
$(l + A $ 1 Z)/2, h = $(1- Aj 1 Z)/2, with A* = $" 1 A$, 
so that if q is invertible one has 



G = 



1 



T + t ' 



(9) 



where G is now seen as the Cayley transform of T = 
A~^Z The fraction symbol is here not ambiguous, 
as we are dealing with commuting matrices, [T — 1, (T + 
l)- 1 ] = [T + 1 - 2 • 1, (T + I)" 1 ] = 0. The matrix 
T = Aj 1 Z can be seen as the orthogonal part of the polar 
decomposition of Z . The latter is the generalization to 
matrices of the polar decomposition of a complex number 
z = |z|e iarg2: . For a generic matrix M one can write 
M = PU, with U unitary and P = VMAft positive 
scmidefinite (note that detM = pe l6 , with p = | detM| 
and e w = det U). If M is normal, [M, Aft] = 0, then one 
can use without ambiguity the notation |Af | = V MA'P — 
VAftAf. In our case, from Eq. (J5J we immediately obtain 
the relation P 2 = ZZ T = A% and thus U = T = $ T *, 
real. Note that, from Eq. ©, the orthogonality T T = 
T _1 readily implies the antisymmetry G T = — G. Since 
G, and thus the ground state \^z), can be expressed in 
terms of T, the latter will be of central importance in our 
analysis. In fact, many of the properties of the ground 
state we are interested in can be related to or derived 
from the properties of T, 

Unless Z is not invertible, T is uniquely defined. Note 
that det 2^0 implies ^ SpA, i.e., the system is gapful. 
In a gapful phase, therefore, T is always well defined and, 
provided T+ 1 is invertible, G exists and the GS is of the 
form J7J. The inverse Cayley transform then gives T — 
(1 + G)/(l - G), which implies detT = 1. Indeed, the 
spectrum of the real antisymmetric matrix G is given by 
complex conjugate pairs of purely imaginary eigenvalues 
(with the exception of an unpaired zero eigenvalue for L 
odd), so that for each eigenvalue 1 + \q of 1 + G there 
exists an eigenvalue 1 — Xq = 1 + Xg of 1 — G. 

The determinant of T, or equivalently the sign of det Z, 
turns out to correspond to number parity, i.e., detT = 1 
and det T = — 1 for ground states of even and odd parity 
respectively. We give a detailed proof in Appendix lAl 
Here it is important to notice that whenever — 1 S SpT 
the matrix T + 1 is not invertible and G cannot be de- 
fined. Equivalently, g is not invertible and Eq. (jSJ) can- 
not be solved. However, with the elementary canonical 
transformations described in Appendix^ one can always 
reduce to a T' such that — 1 ^ SpT', so that the ground 
state can be expressed by Eq. in terms of the trans- 
formed quantities. This will be useful in deriving the 
fidelity formula of Subsec. TlIBI In addition, provided L 
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is even and detT = 1, one can write the ground state in 
a form which is valid independently of the presence of —1 
in SpT and can hence be considered more general than 
Eq. 0. This is 



^4i[cos(^/2)|00),,_, + sin(^/2)|ll) 

v, — v\ 



(10) 

where {e ±l6 "} L J= x = SpT and \0) u denotes the vac- 

uum (occupied) state of the new pairs of fermionic modes 
c v = c' 2u _ l , c^ v = c' 2 „, with c'j = Rf k Ck- The matrix R 
is the one putting T and G in block diagonal form (see 
Appendix EJ. In Ref. @ Eq. JTDJ was derived starting 
from Eq. Q) , but it should be clear from the discussion in 
Appendix ^ how to prove it for cases where G cannot be 
defined. One also recognizes when the Ansatz (7J is not 
valid: if — 1 S SpT then 9 V — ir for some v — 1, . . . , L/2 
and hence the corresponding term in the direct prod- 
uct (|10|l is simply lll)^^^, which can be described by 
expression (J7J only in a limiting sense (the weight of the 
jOO)^-!/ term cannot vanish unless the normalization J\f 
diverges) . 



B. Ground state fidelity 

We are interested in characterizing how the ground 
state changes when the set of parameters of Eq. are 
slightly modified, i.e., when (A, B) -> (A + 8 A, B + SB) 
or equivalently when Z — > Z = Z + 5Z . In order to 
characterize these changes we use the the ground state 
fidelity, which was introduced in Ref. and discussed 
in Ref. 0: 



^(Z,Z) = |<*z|*z; 



(11) 



As already pointed out in Ref. [8j , in order to evaluate it 
we can resort to the theory of coherent states ^jj. In fact, 
the GS construction given by Eq. (J7J) coincides with the 
construction of a coherent state for the group SO(2L, H). 
In general the construction of a coherent state relative to 
a Lie group Q involves the action of an element of a uni- 
tary representation of the group over the so called max- 
imal weight vector of the Lie algebra A. In our case, we 
have that the group coincides with SO(2L,Wi) and the 
relative Lie algebra so(2T,IR) is isomorphic to the alge- 
bra generated by the operators {c{Cj, cjcj, i(cjcj — cjcj)}. 
The group SO(2L, H) has two unitary irreducible repre- 
sentations over the 2 L fermionic Hilbert space that splits, 
according to these representations into two irreducible 
subspaces even and odd, i.e., H.f = Wj(J)7i^, hav- 
ing the basis vectors with an even and odd number of 
fermions respectively. If we now use the even representa- 
tion, we can generate the generic coherent state by acting 
with an element U + = cxp(A), A € so(2T,IR) upon the 
maximal weight vector of the algebra that in this case is 
the already introduced fermionic vacuum (cj|0) =0). If 
we now use (A)^ = —GijCjc] we have that the generic 



coherent state can be identified by a skew-symmetric ma- 
trix G (in our case determined by Z) and has exactly the 
form Q. This identification allows us to use the general 
result for the scalar product (overlap) between fermionic 
coherent states defined by the antisymmetric matrices G 
and G, namely 0] 



det(l + G*G) 



1/2 



det(l + GtG)!/4 det(l + GtG) 1 / 4 



(12) 



We can then explicitly write the ground state fidelity in 
terms of T. In fact, as described in Ref. Q one finds 



T(Z, Z) 



\ 



det 



T + T 



1 +T~ 1 T 



I det 



(13) 



It is not difficult to get the latter result by substituting 
G = (T - 1)/(T + 1) into Eq. $2&. Whenever one of 
the ground states cannot be expressed directly into the 
form 0, two possibilities are present. Either the ground 
states belong to different sectors of fermion number par- 
ity (that is detT = -detf and de^T^f) = -1) and 
the fidelity vanishes, as correctly reproduced by Eq. Qliip. 
or they belong to the same sector but — 1 is an eigenvalue 
of T and/or T. In the latter suitable canonical 

transformation can always be found (similar to those ex- 
plained in Appendix lA"|l such that (i) —1 is absent from 
the spectrum of both the transformed matrices T' and T' 
and (ii) det(T' + f') = \ det(T + f )|. Consequently, one 



has T = ydet[(T' + T')/2] = ^\ det[(T + T)/2]| and 
formula (|13|l is proven for all the possible cases. 

We also recall that if T~ lr f E SO L {WC) and L is even 
then 



L/2 



^(Z,Z) = JI|cos(e„/2)| , 

1^=1 

where SpC^T ) = {e ±lQ -}^i\. 



(14) 



C. Fidelity derivatives 

The fidelity function compares ground states at dif- 
ferent points in the parameter space, giving a measure 
of their 'distance'. Intuitively, in a critical region the 
ground state changes significantly even for small varia- 
tions of the Hamiltonian parameters. In other words, 
the rate of orthogonalization diverges in the proximity of 
QPTs. It is then reasonable to characterize the degree of 
criticality by the rapidity of the ground state variation 
as a function of the system parameters. In practice, this 
is clearly related to derivatives of the fidelity function. 

We fix a point Z in the parameter space, with the 
corresponding ground state \^z)- The fidelity between 
\^>z) and the GS l^^)' corresponding to a second point 
Z, can be regarded as a function of X = Z — Z, defined 
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by F Z (X) = T{Z,Z). Clearly, the function F Z {X) has 
a maximum for X = 0, where F z (0) — 1. Consequently, 
the first derivatives of F Z (X) with respect to X must 
vanish. In order to analyze the fidelity behavior around 
X = one has then to consider second derivatives, i.e., 
the Hessian matrix of F Z (X) calculated at X — 0, which 
we denote by Hp z (0) . The eigenvalues of the Hessian ma- 
trix correspond to the principal curvatures of the hyper- 
surface generated by the function F Z (X) at X = 0. The 
eigenvector corresponding to the lowest (highest in mod- 
ulus) eigenvalue gives the direction of most rapid change 
in the GS. In the proximity of a QPT one expects a (neg- 
ative) divergence of the smallest Hessian eigenvalue. We 
then use h(Z) = min{Sp[i?F z (0)]} as a measure of the de- 
gree of criticality of the point Z . For the two-dimensional 
parameter space /i-7 considered in this paper one has 
Z = Z(ji,i), X = X(5n,6"/), and 



h(Z) 



[dLFz 



d 2 Sl F z - 



(15) 



(dLFz 



dLFzY 



±{d Sll d 5l F z y 



where all the derivatives are calculated at 8)1 = 5"/ 
(corresponding to X — 0). 



III. VARIABLE RANGE FERMIONIC XY 
MODELS 

We begin this section by briefly reviewing some aspects 
of the traditional XY spin chain in a transverse mag- 
netic field 0|- I n accordance with recent literature, we 
write the corresponding ferromagnetic Hamiltonian as 



H = 



E 

3 



1+7 



1-7 



V V 



where crj* for a 



(16) 

x,y,z are the usual Pauli operators 
relative to the j'-th site, 7 is the anisotropy parameter, 
and A defines the strength of the magnetic field in the z 
direction. Actually, in the original papers 0, IT^j an an- 
tiferromagnetic spin chain was chosen. As recalled later, 
however, this makes no difference for the ground state 
phase diagram. 

In Eq. (|16fl we did not specify the range in the summa- 
tion signs. Two possible choices are indeed possible, (i) 
the cyclic chain with periodic boundary conditions and 
(ii) the free-ends chain. We will return on this point when 
considering the fermionic version of this Hamiltonian. 

In fact, by using the Jordan- Wigner transformation, 
Eq. I|16l) can be rewritten in terms of the fermionic oper- 
ators 



(17) 



where 



<7j ± «<rJ)/2 are the Pauli raising (+) and 
lowering (— ) operators. For a free-ends chain with L sites 



one gets 



L-l 



H = 



L 

]{c]c J+1 + 7 c]c] +1 + h.c.) - J2 *(2ct Ci - 1) 



/_^V-J-J-|-J. ' i-J-J-t-l ■ ' 

i=i 3=1 



while for the cyclic case additional boundary terms ap- 
pear, which are however negligible in the thermodynamic 
limit 0- Clearly, apart from a constant term, Eq. 118|) 
is a special case of Eq. (QJ. 

The Hamiltonian (|18fl corresponds to a free-ends chain 
in terms of the Ci fermionic operators. Analogously, 
one can consider the "c-cyclic" problem, where periodic 
boundary conditions are applied to the fermionic chain. 
In Ref. |9j the c-cyclic problem is shown to be equivalent 
to the periodic spin chain in the thermodynamic limit. 
Here, we generalize the fermionic XY Hamiltonian by 
introducing variable range couplings and studying both 
the c-cyclic and the c-free-ends problems. In the free- 
ends case, the Hamiltonian with a coupling of range r 
reads 



H 



L — X niin(r,L 

E E 

j=i *;=i 



-3) 



(4c 



-h.c. 



L 

3=1 



(19) 

where [i plays the role of the chemical potential. For 
r = 1 one reduces to the nearest neighbor coupling of 
the XY model (apart from a constant, the identification 
with Eq. H18|) takes place by putting \i — 2A), while for 
r = L — 1 one gets the fully coordinated fermionic system. 
It is worth noticing that the spin Hamiltonian obtained 
by the inverse Jordan- Wigner transformation of Eq. 
is not the variable range extension of the XY model 
For example, the fully coordinated spin system, known 
as the Lipkin-Meshkov-Glick model [lg , does not corre- 
spond to the fully coordinated fermionic graph discussed 
below. We also point out that, apart from a constant, 
the sign of the Hamiltonian is changed by the canonical 
transformation c'j = ct, so that the GS phase diagram 
of —H is equal to that of H . This is why both the fer- 
romagnetic and antiferromagnetic coupling in the XY- 
model give rise to the same phase diagram (except for 
the sign of the magnetic field) . 

It is easy to express our variable range Hamiltonian 
in the form of Eq. QJ. Since in the Hamiltonian the 
couplings only depend on the distance between sites, A 
and B are Toeplitz matrices, i.e., the element of the j-th 
row and the fc-th column depends only on the difference 
j — k. In the free-ends case, by using the discrete step 
function 9 one can write 

A jk (r) =-[(!*- + 0(r ~ \j - k\)] , (20) 
Bjk(r) = - 7 sign(fc - j)9(r - \j - fc|) , (21) 

where j,k = 1, . . . , L, r = 0, . . . , L — 1, and we adopted 
the usual conventions 9(0) — 1, sign(0) = 0. In the case 
of periodic boundary conditions, where translational in- 
variance modulo L takes place, A and B become circu- 
lant, i.e., each row is a circular shift of the previous one. 
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For example A21 = Ail, corresponding to the fact that 
the weight of the term c^Cj in the Hamiltonian must be 
equal to that of the term c 
can then be written as 



1 U L 



The matrices 



B jk {r) 



-[^-l)5 jk + 6(r-\j-k\) + 
-6(\j-k\-L + r)} , (22) 
-7sign(fc -j)[0(r - \j - + 
-e(\j-k\-L + r)} , (23) 



where now r — 0, . . . , [(L — 1) /2J , with [• • ■ J denoting 
the 'floor' of a real number. The cyclic case r = L/2 for L 
even (corresponding to the fully coordinated system) has 
to be treated separately: in this case Aj k = — [(/i— l)<5ufc+ 
1, while B jk = B ]k {r = L/2 - 1), i.e., B hL/2+1 =0 0- 
In order to visualize at a glance the simple structure of 
the above matrices, we write explicitly the matrix B for 
the case r = 2 with periodic boundary conditions, which 
is 



B 



1 1 

1 '•• 







-1 



-l\ 



- 1 





V 





1 '•• 
1 1 







(24) 



0-1-1 / 



We discuss in some detail the cyclic case. Circulant 
matrices have some useful properties which make their 
analysis very simple: (i) they all can be diagonalized by 
the same unitary transformation and (ii) their eigenval- 
ues are given by the discrete Fourier transform of their 
first row. Property (i) implies that circulant matrices 
commute and that sums and products of circulant ma- 
trices are also circulant (as the transpose as well). Fur- 
thermore, whenever a circulant matrix is invertible, the 
inverse is also circulant. In brief, circulant matrices form 
a commutative ring. The diagonalizing unitary matrix is 
usually written as 



, (circ) I 



■»2ir0'-l)(fc-l)/i 



(25) 



where j, k — l,...,L. The eigenvalues, according to 
property (ii), are instead given by 



A 



(circ) _ 



E 

k=l 



m lk e 



-<27r(j-l)(fe-l)/i 



(26) 



where mj k are the entries of the circulant matrix and 



again j,k = l,...,L. Note that A^ clrc ' ) = YlX=i' mik ^ 



while, in the case of real matrices, X^"f + 2 = .\( C1IC )* for 
j = 2, . . . , L, where the star denotes complex conjuga- 
tion. Consequently, whenever the circulant matrix is real, 
one has \_(L — 1)/2J pairs of complex conjugate eigenval- 
ues, with the additional unpaired real eigenvalue A^ " ^ 
for L odd and the two additional unpaired real eigenval- 
ues \^i" c \ \^J^ +l for L even. For symmetric matrices, 
where eigenvalues are real, complex conjugate pairs be- 
come degenerate, so that the columns of matrix Ij25(l can 
be linearly combined in pairs to get real eigenvectors. 

The diagonalization of A, B, Z = A — B is then 
straightforward. We denote by a, £ IR, /?, £ ilR, 







Pj the eigenvalues of A, B, Z respectively 



(j = 1, . . . , L). Since Z is real, < £ Sp(Z) => C* e Sp(Z). 
Clearly, in the polar decomposition Z = A$T the posi- 
tive part A$ = \J ZZ T = \Z\ and the unitary part T are 
also circulant. If Z is invertible, T = Zj\Z\^ x and its 
eigenvalues are simply given by tj — e l6j — Cj/\Gj\> i- e -j 
9j = arg(Cj) = arccos(3?Cj/IOI)- Then from Eq. {Bj one 
gets 



1 + e i03-9j) 



M+l 

n 

J=2 



(27) 

where Qj = 8j — 8j, M = [(L — 1)/2J is the number of 
complex conjugate pairs, and Sq is a function which is 
zero or one depending on whether @i (and Ol/2+i for 
L even) is it or zero (modulo 2-7r) respectively. Hence, 
either all of the unpaired L — 2M eigenvalues of T~ X T 
are equal to one, or the fidelity vanishes, corresponding 
to a parity change in the ground state. 

We evaluate explicitly the eigenvalues for the variable 
range c-cyclic case. By using Eq. I|26|) one finds 



Cj(r) = M+(-!)^r,L(i+i)/2j + 

2TTk(j - 1) 



(28) 



k=l 



L 



. . 27rfc(j - 1) 
ij sin 



where r = 1, \(L — l)/2], with [~. ..] denoting the 
'ceiling' of a real number. The function 8 r< 1 (l+i) /2J i s 
non-zero only for r = L/2, with L even. The sums in 
Eq. H28|) can be usefully rewritten as Y^k=i cos(k9) = 
sin(r6>/2)cos[(l -I- r)0/2]/sin(0/2) and ELi sin ( fc61 ) = 
sin(r-0/2)sin[(l + T-)0/2]/sin(0/2). The eigenvalues of the 
matrix Z then assume the following expression 

-Cj(r) = M+ (-l) J 5 r ,L(L+i)/2j + 



(29) 



} sm(rKj/2) 
sin(Kj/2) 



where Kj — 2tt(J — 1)/L and again j — 1, . . . , L. For 
the nearest neighbor range r = 1 one recovers the well- 
known result of the XY-model, which, in terms of the 



single particle spectrum, yields Aj(r 
2^j (cos Kj + fi/2) 2 + 7 2 sin 2 Kj $22} . 
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IV. FULLY COORDINATED FERMIONIC 
HAMILTONIAN 



In this section we discuss the long range fermionic sys- 
tem, where, apart from the diagonal elements, A is the 
adjacency matrix of the complete graph. The entries of 
the matrix B can instead be interpreted as the orien- 
tations of the couplings, so that we have the structure 
of a directed graph. To be consistent with the nota- 
tion of Ref . in the following we change the sign of the 
Hamiltonian with respect to the previous section, so that 
A — > —A, B — ► —B, Z — > —Z. As noticed above, this 
does not affect the phase diagram. 



6 fx, A2 = £7 and one has 



d^Fz 



dsuds-yFz = 



S , 



_2_ 

4 

0»-i) 2 
4 

(M-l)7 



s 



s . 



with s = T,f=2 ls J and s j = (^0)7101 

quently, one straightforwardly finds 



ft[^(/i,7)] 



1 



[(M-1) 2 +7 2 ]S, 



(34a) 
(34b) 
(34c) 

Conse- 
(35) 



A. Cyclic case 

We proceed by considering the cyclic fully coordinated 
case r = \(L — l)/2] . One then finds £1 = fi + L—1 and 



while the other Hessian eigenvalue is always zero and 
hence det Hp z (0) = 0. We also note that, apart from the 
discussion on first order QPTs, all these results apply 
equally well to the even L case. 

We are interested in studying h in the limit L — > 00, 
where a divergence of h reflects a fast fidelity drop. We 
£. = fj, — l 4- i<y[i + (— 1) J ] cot — (L even) , (30) therefore analyze the behavior of the functions sj enter- 
ing the sum S. For (1 ^ 1,7^0 none of these functions 
is singular, independently of L, so that in the thermo- 
dynamic limit one can approximate their sum with an 
integral 



(1 — 1 + «7 



2 4 



i + (-iy +1 



tan ■ 



(L odd) 



(31) 



for j — 2, . . . , L, where Kj = 2ir(j — 1)/L. Hence, con- 
cerning the imaginary part = i(3j, for L even one 
has 3C2fe = 27Cot[7r(2fc — 1)/L] and 3C2fc-i = 0, while 
for L odd SJ&fc = 7Cot[7r(2fc — l)/2£] and ^C,2k-i = 
— 7tan[7r(2fc — 2)/2L\. For the real part — aj, we 
note that the above simple result can also be readily ob- 
tained by considering the expression 



A = L|^ 1 )(^| + (/i-l)l = 
= {L + (i -l)|Vi>(Vi| + (/*■ 



1)(1 



(32) 
l^iX^il) , 



where |^>i) is the uniform vector given by the first column 
of the matrix ip cm defined in Eq. I|25|l . in components 
^jcnc) _ Yjyjjj, Equation l|32|l directly diagonalizes A, 
i.e., SpA = {L + (i — — 1}, where the eigenvalue 
a — (i — I is L — 1 times degenerate. 

To proceed with our analysis we now fix L odd, as 
usual for the XY-model. We express the fidelity in the 
product form of Eq. (|27J) and calculate the derivatives 
entering Eq. I|15|l . We neglect the function 5® in Eq. (|27|l , 
as it only gives rise to the first order QPT placed at 
(1 = 1 — L, which is uninteresting in the thermodynamic 
limit L — ► 00. Then 



d 2 F z 



dXhdXi 



1 ^ 96, 



\ h =o 



4 ^ 

J=2 



\ h =o 



d\i 



(33) 



Ah=0 



where the parameters {Xh}h=i,2 are the variations Ai 



1,7^0) 



tan 2 x 



l r /4 

Wo t[Gu-l) 2 + (7tan.x)212 



cot 2 X 



[(/i-l) 2 + ( 7 cota;) 2 ] 2 
1 L 
4| M -1|| 7 |(| M -1| + | 7 |) 



da; 



(36) 



from which it follows the divergence h oc L. This corre- 
sponds to the fact that, for any value of the Hamiltonian 
parameters, the fidelity between different ground states 
vanishes in the thermodynamic limit. 

A different situation takes place for fi = l,j ^ 0. In 
this case one has s 2 j+i = cot 2 (nj/L)/^ ~ L 2 /(nj-f 2 ) 2 
for L —> 00, so that it is readily seen that the sum S 
scales at least as L 2 . For large L we find 



, O(L) 

1 3 = 1 

1 ° (L) / T \ 2 \J 2 
1 V^J/ 6 7 



(37) 



were in the last summation we used O(L) — > 00 and hence 
the well known result for the Ricmann £ function £(2) = 
Yl^jLiJ 2 = tt 2 /6. One then finds that the line (i = 1 
corresponds to a faster decay of the fidelity in the limit 
L — > 00. This puts in evidence the presence of a second 
order QPT, similarly to the case of the XY-model Q for 
(i = 2 (i.e., A = 1). The same analysis can be repeated 
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FIG. 1: Fidelity :F min = mm[T{Z, Z Sfl ), T{Z, Z 5l )] (see 
Ref. 0) in the /i-7 plane for L = 1001, with Sfi — — 0.1. 
At the critical lines fi — 1 and 7 = one observes an evident 
drop. 




FIG. 2: Decimal logarithm of the absolute value of h(Z) in 
the /i-7 plane for L — 1001. Note the singular behavior at 
fi = 1,7= 0, where the function is non-analytic. The critical 
lines fi — 1 and 7 = emerge in a clear way. 



for m * 1,7= 0, where we get 5 - (l/2)[£/(/i - l) 2 ] 2 . 
Summarizing, one has 



(M-l) 2 +7 2 



Mm = 1,7*0) - 



16| M -l|| 7 |(|M-l| + |7l) 2 
2 



1,7-0) 



1 / £ 



/i — 1 



, (38a) 
(38b) 
(38c) 



Finally, we note that the point ft — 1,7 = has a pecu- 
liar singular behavior. There, all but one (£1 = L + fj,— 1) 
the eigenvalues of Z vanish and the unitary T is unde- 
fined (as at a first order QPT where level crossing takes 
place). Consequently, the angles 0» cannot be evaluated 
and the above formalism does not apply. 



The analysis of second order QPTs can independently 
be done in terms of the derivatives of the GS energy E — 
Tr(A- A)/2 = Ei=i(^C? ~ IGD/2. Except on the fx axis, 
where Eq — for fi > 1 and E = [L— l)(/x— 1) for /i < 1, 
the energy scales in a superextensive way, due to the fully 
connected nature of the system. Indeed, while TrA = 
LfJL, for 7 ^ one has Tt\B\ < TrA < Tr|A| + Tt\B\ ~ 
Tr|_B| cx L\vlL for L — ► 00. This makes it impossible 
to define the energy density Eq/L in the thermodynamic 
limit. The density of energy derivatives can however be 
calculated for finite L and in the TDL as well. All the 
second derivatives are finite and continuous for /1 ^ 1, 
7^0. For fi 7^ 1, 7 = the second derivatives d^d^E^/L 
and d^Eo/L are analytical for L finite, but in the TDL 
one finds the discontinuous behavior <9 7 d p i?o(/i ^ 1,7 — > 
O^/Z - ±(1/tt)/(/x - 1) and the divergence 9 2 S (^ * 
1,7 = 0)/£ cx S(jt ^ 1,7 = 0)/L oc L It is hence 
reasonable to identify the critical line 7 = with a second 
order QPT. Note, however, that the gap is finite at the 
critical point. The discontinuity and the divergence of 
the second derivatives of energy and fidelity at this gapful 
transition in the TDL is due to the divergence cx L of 
some of the single particle energies Aj — This is a 
consequence of the unphysical long range nature of the 
fully coordinated model. Indeed, from Eq. (|26|l one has 

101 < ELi l^i*l < ELi(IAfe|+|£i fe |), so that 101 cx L 
only if r cx L. Another second order QPT is found at 
fi = 1, 7 7^ 0, where d^Eo/L cx \nL, while the other 
second derivatives are continuous. Here, however, one 
also has the vanishing of the lowest single particle energy, 
since $tQ = for j = 2, . . . , L and 3(fc — > for L — > 00 
with k odd. 

In Fig. n we report the GS fidelity in the /i-7 plane, by 
using JF min = min[.F(.Z, Zsp), F{Z, Zg-y)], where Zg^ = 
Z{p + Sfi, 7) and Zg 1 = Z(/j,, 7 + £7), as in Ref. 8]. This 
corresponds to plotting the minimum of the fidelity with 
respect to variations along the two axes, which are related 
to dg Fz and <9? Fz- Since the critical lines are parallel 



to the axes, these variations are sufficient to fully char- 
acterize the phase diagram. A detailed analysis of all the 
second derivatives, including also ds^ds^Fz, is instead 
provided by the study of h(Z). The behavior of h{Z) 
in the /i-7 plane is shown in Fig. [21 for L = 1001. The 
qualitative agreement with Fig. ^ is evident. 

Equations (|38|l give the scaling of h in the thermody- 
namic limit. It is also interesting to look at the finite-size 
scaling of this function, in connection with the next sec- 
tion about the free-ends fully coordinated system, where 
the simple analytical results obtained above for L — > 00 
are not available. The finite-size scaling for the critical 
line fi = 1 is shown in the upper panel of Fig. [3J The 
peak of h at fi = 1 increases (in modulus) with L 2 , as 
expected from Eq. l|38b|) . Away from the minimum, the 
function is instead proportional to L, as from Eq. (|38a|l . 
Consequently, the different scaling behavior makes the 
peak shape progressively narrower as L increases. It is 
possible to plot the function h in rescaled units in order 
to put in evidence these qualitative considerations. By 
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-0.02 > ' 1 ' 1 

-10 10 

(M)L 



FIG. 3: Finite-size scaling (upper panel) and data collapsing 
(lower panel) at the critical line fi = 1. The upper figure 
shows the enhancement of the negative peak of h at fi = 1 
for 7 = 1.5 and L = 101, 201, 1001. In the inset the 
curves are plotted after dividing ft by L, which allows a direct 
comparison with the asymptotic curve (dashed line) given by 
Eq. (I38al . The lower figure replots the same curves of the 
upper figure in rescaled units. Most of them are practically 
indistinguishable and give rise to the thick line in the plot. 



dividing the peak depth by L 2 and rescaling the distance 
fj, — 1 from the critical point by L (therefore compen- 
sating the peak narrowing) one gets a curve practically 
independent of L (data collapsing) , as found in the lower 
panel of Fig. [3] This can be understood by looking at 
Eq. i|38a[) in the proximity of the critical line /i = 1, i.e., 
for \fi — 1| -C 7, which is well verified for the parameters 
of Fig. Then one finds h ~ — L/16\fx — so that, 

for a given value of 7, the quantity h/L 2 turns out to 
depend only on the combination \fx— 1\L. Similar results 
were observed for the critical line 7 = 0. 

We conclude the analysis of the cyclic fully coordi- 
nated system with a couple of additional remarks. As 
anticipated above, det Hp z (0) = 0, as one eigenvalue 
identically vanishes. The components of the eigenvector 
v = (v^, v~f) associated with this zero eigenvalue are read- 
ily found to be related by w 7 = -{dj F z /d Sf _ l d s .yF z )v ll . 



Hence, since -dj F z / d Sf _ l d s ~ ( F z = j/(fi- 1), the eigen- 
vector is directed in the radial direction from the point 
fj, = 1,7 = 0. By recalling that the two eigenvectors of 
the symmetric matrix Hp z (0) are mutually orthogonal, 
one immediately concludes that the direction of fastest 
GS variation is the azimuthal one with respect to this 
shifted origin. Finally, we comment on first order QPTs 
for L even. In this case one has two first order criti- 
cal lines, namely /j, = 1 — L and fi = 1, corresponding 
to the vanishing of the eigenvalues £1 — L + /1 — 1 and 
Cl/2+i =fi-loiZ. 



B. Free-ends case 

In the following we discuss the free-ends fcrmionic 
graph. Some results on this system were already re- 
ported in Ref. @. In this case the relevant matrices are 
not circulant and the analytic results used before are not 
available. In particular, it is worth stressing that the 
matrix Z in the free-ends case is in general not normal, 
i.e., [Z, Z T ] 7^ 0. We therefore proceed by calculating the 
phase diagram numerically. 

Let us first recall some analytical facts ||. For 7 = 
the resulting number-conserving single-particle Hamilto- 
nian is the same as for the cyclic case. For (/j, = 0, 7 = 1) 
the matrices Z and Z* = Z T become instead lower and 
upper triangular respectively. By explicit computation 
one finds (ZZ T \j — 4min(L — i, L — j), which has the 
last column and row identically vanishing. Accordingly 

S Sp \Z(0, 1)| V L. H3|. We also recall that changing 
the sign of 7 simply corresponds to transforming Z into 
Z T . Hence, since Sp(ZZ T ) = Sp(Z T Z), the matrix A is 
not affected by the transformation. Then one can define 
\?(— 7) — $(7) and consequently $(— 7) = ^(7), where 
the matrices $ and ^ were introduced in Sec. [HJ This 
implies that T(- 7 ) = $ T (-7)*(- 7 ) = [$ T (7)*( 7 )] T = 
T T (7), so that the overlap behavior is symmetric with 
respect to the 7 = axis. 

We start from the analysis of first order QPTs. Dif- 
ferently from the cyclic case, the position of these tran- 
sitions in the parameter space depends on L in a non- 
trivial way, reaching however a simple configuration in 
the thermodynamic limit 8] . First order QPTs are given 
by level crossings which take place when the lowest sin- 
gle particle energy (the gap) exactly vanishes, i.e., when 
det A = I det Z \ = 0. Typically, at these points one also 
has a discontinuous change in det T, which jumps from 

1 to —1 or vice versa (corresponding to a sign change in 

detZ), implying T = det[{l + T- 1 f)/2] = when- 
ever T and T are calculated at opposite sides of the tran- 
sition. Strictly speaking, at the transition point where 
det Z = the unitary T is instead undefined, similarly 
to the azimuthal angle of polar coordinates at the origin. 

In order to calculate the position of first order QPTs in 
the parameter space one then has to solve the equation 
det Z = 0. The corresponding curve in the /z-7 plane 
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-2-1012 -2-1012 
li l» 

FIG. 4: Left panel: zeros of det Z in the fi-j for L = 3, 
given by 7 2 = — (/i — l) 2 (/i + 2)/(3/i — 2). Right panel: phase 
diagram in the /i-7 plane in the limit L — > 00 for L odd. Solid 
and dashed lines correspond to first and second order QPTs 
respectively. 



depends on L. Furthermore, its thermodynamic limit 
depends on the parity of L. 

The results for L even were already reported in Ref . || . 
For L = 2 the model is trivially solvable and one finds 
det Z = /1 2 + 7 2 — 1, so that the critical line is the circle 
of radius one centered at the origin. For L = 4 and L = 6 
the equation det Z — is of second and third degree in 7 2 , 
respectively, so that relatively simple analytical solutions 
are available. One gets in this way a sequence of closed 
curves crossing the n axis at fi = 1 — L and fj, = 1, while 
the 7 axis at 7 = ±1. Numerically, in the limit L — > 00 
one finds the boundary 7 = ±1, /i < 1 and I7I < 1, /i = 1. 

For L odd the results are different. The equation 
det Z = for the case L = 3 is again readily solved ana- 
lytically. For any L, the resulting curve is now open and 
made of two branches. One branch starts at ji = 1 — L, 
passes through ji = 0,7 = ±1, and diverges (as a func- 
tion of n) at n — fi c , where /i c — 2/3 for L = 3 (see left 
panel of Fig. ■ The second branch trivially reduces to 
the point (J, = 1,7 = 0. The numerical analysis yields 
fi c — > 1 in the thermodynamic limit, so that the critical 
boundary (apart from the singular point /j, = 1,7 = 0) 
is given by 7 = ±l,/x < 1 and I7I > 1,/i = 1 (see right 
panel of Fig. 0}. 

The fact that the asymptotic first order boundaries are 
different in the even and odd L cases is at first sight sur- 
prising. In the limit L — ► 00, where L + 1 ~ L, one 
would expect to find the same result in the two cases. 
However, the phases separated by the considered QPTs 
just differ by the fermion number parity. The correspond- 
ing level crossings exchange energy states with different 
parity (recall that the number parity is conserved by the 
Hamiltonian) . Loosely speaking, these states, although 
orthogonal for any L, become less and less different as L, 
and hence the average number of fcrmions contained in 
the GS, increases. In conclusion, in the thermodynamic 
limit such first order QPTs are of little physical inter- 
est, in the sense that in practice it would be very difficult 
to discriminate between such phases, as it would require 
the capability to distinguish between L and L + 1 for 
L> 1. On the other hand, the different asymptotic be- 




FIG. 5: Finite-size scaling (upper panel) and data collapsing 
(lower panel) at the critical line /1 = 1 for 7 > 1 and L 
even, where (only) a second order QPT is present. The upper 
figure shows the enhancement of the peak of h at /1 = 1 for 
7 = 1.5 and L = 100, 200, . . . , 1000. In the insets, the law 
1 — A'min ~ 1/L for the peak position (error bars correspond to 
numerical precision) and the function h/L are shown. In the 
lower figure,the same curves of the upper figure are replotted 
in rescaled units, practically merging into an indistinguishable 
thick line as in the cyclic case. 



havior depending on the parity of L is not unreasonable 
for a quantity which is in turn related to a parity, namely 
the fermion number parity. 

We complete the analysis of the phase diagram by look- 
ing at the fidelity. First order QPTs appear as sudden 
drops of the fidelity, which falls to zero when the consid- 
ered GSs are in different phases. Besides these discon- 
tinuous drops, one also observes a smooth but evident 
fidelity decay in the proximity of the lines /j, — 1 and 
7 = 0, which we identify with second order transitions 
|g. Hence, apart from first order QPTs, the phase di- 
agram of the free-ends case corresponds to that of the 
cyclic case. 

It is worth discussing in some detail the critical line 
fi = 1. As explained above, for L odd and \j\ > 1 a 
first order QPT is also present. For a given value 70 with 
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I70I S (l,oo) and for finite L the (first order) transition 
point takes place at a critical /J, c ("fo) < 1 which tends 
to 1 in the thermodynamic limit (see Fig. 0J. At the 
point /x c (7o),7o the fidelity function is discontinuous and 
its derivatives, and hence the Hessian matrix, are not 
defined. However, everywhere else on the line 7 = 70 the 
eigenvalue h(Z) = mm{Sp[Hp (0)]} can be calculated. 
This allows to analyze the finite-size scaling behavior of 
h in the neighborhood of /z = 1 even for I7I > 1. One 
then finds the same divergence oc L 2 as for second order 
QPTs in the absence of first order transitions. Only, the 
constant of proportionality is different in the two phases 
(i.e., for even or odd number of fermions). For I7I < 1, 
instead, one only has the second order transition and h 
is symmetric around the critical point. The position of 
the minimum of h in the thermodynamic limit gives the 
second order critical line. For finite L this position is 
given by the line /x = /i m i n , independently of 7. In the 
thermodynamic limit the numerical analysis gives 1 — 
/Zmin ~ so that fi = 1 is the asymptotic critical line. 

For L even, the first order QPT interesting the /i = 1 
transition takes place for |7| < 1. Except for this differ- 
ence, everything goes as in the odd L case. The finite- 
size scaling features of h around fj, = /x m ; n for L even and 
7 = 1.5 are summarized in Fig. [S] It is also worth notic- 
ing that, albeit only in the TDL, the gap also vanishes 
in the part of the /i = 1 critical line where the first order 
QPT is absent (e.g., for I7I > 1 with L even). 

Another second order QPT is found along the line 7 = 
0. Here the gap is always finite, apart from the point 
/i = 1 (see the cyclic case). Nevertheless, the finite-size 
scaling of h puts in evidence the critical behavior. The 
position of the minimum of h is here independent of L, 
being always given by the line 7 = 0. 

Some insight in the nature of the latter transition 
can be gained by looking at the perturbative expan- 
sion of the ground state around 7 = 0. First of all 
we use a series expansion in 7 to prove analytically the 
L 2 divergence of h in 7 = 0. Along the line 7 = 0, 
apart from the first order QPTs located at // = 1 — L 
and n = 1, the ground state is constant. Therefore 
d 5^i F z = dsfidg-fFz = and consequently h = d^F z . 
In practice, one can reduce to the one-dimensional pa- 
rameter space given by 7 and use the series expan- 
sion in terms of ordinary derivatives presented in y|, 
T[Z{\),Z{\ + SX)] = 1 + (l/16)Tr(iT(5A) 2 + 0(<5A 3 ), 
where we assumed det T = 1 so that T = e K and 
the prime denotes derivation with respect to the one- 
dimensional parameter A (not to be confused with the 
magnetic field strength in the XY-model). For simplicity 
we consider the case fj, > 1, where T(7 = 0) = 1 and con- 
sequently T'( 7 = 0) = T( 7 = Q)K'(>y = 0) = #'(7 = 0). 
Lengthy but simple calculations (see Appendix 0) yield 

+ 2L(2n - 3) + 4(/i - l)(3/x - 5)] , 

(39) 



so that one finds h(-y = 0) = (l/8)Tr[T'( 7 = 0)] 2 ~ 
— (l/24)[i/(/i — l)] 2 in the thermodynamic limit. Again 
by series expansion one has G(-f) — T'(-f = 0)7/2 + 
0(7 2 ). The divergence of the elements of T" at 7 = 
(see Appendix [BJ) in the thermodynamic limit therefore 
accounts for the rapid orthogonalization rate around the 
critical line given by the [i axis. 



V. SINGLE SITE ENTANGLEMENT 

The idea of systematically using measures of quantum 
correlations in order to characterize QPTs has been ex- 
plored in the past years for general models \A and in 
particular for the case of fermion systems 0, . It was 
found that some of the proposed measures are able to 
spot, with their peculiar behavior, the undergoing QPTs; 
in particular, in some cases the behavior of their deriva- 
tives is universal, in the sense that they diverge with the 
correct critical exponents. Here we would like to com- 
plete the fidelity analysis of QPTs in the cyclic fully con- 
nected system of Subsec flV Al bv studying of one of the 
simplest of such measures: the single site entanglement 
Si = -Tr(p i log 2 p i ), where Pl = Efc,fee{o,i} Ph,k\h) (k\ 
is the 2x2 single site reduced density matrix [l^. It 
can be easily shown that for the system under study 
we have that po,o = 1 — Pi.i = where n is the 
density. Thus the single site entanglement is given by 
Si = —n log n — (1 — n)log(l — n), where here the log- 
arithms are taken in base 2. It can be explicitly com- 
puted by noticing that n = d(Eo/L)/dfi. As already 
described, in our case we have Eq = (TrA — TrA) /2 = 
(Lfi — J^j IOD/2 and n can be obtained by calculating 
the TDL of \(j\)/L = T u . The last equality fol- 

lows from the following observations: (i) since T is cir- 
culant all the elements of its diagonal are equal and each 
of them corresponds to Tu = (TrT)/L = Tj/L, where 
we recall that Tj — Cj/\Cj\ are the eigenvalues of T; (ii) 
if Q € Sp(Z) then Q € Sp(Z) and thus T u £ M; (iii) 
the explicit formulas for Q given in Eqs. (|30|l - l|31|) allow 
to write t\ = 1, for fi > 1 — L, which is always true 
in the TDL, and d^Q] = (/x - 1)/|0| - 3?0/lCj| for 
j = 2, . . . , L. We can thus write 

1 1 L r) TrA 

^ = Z + zE^/IOHl-2^-— . (40) 

The above derivation of Si in terms of Tu shows the phys- 
ical link between the involved quantities and in particu- 
lar with the single particle eigenvalues Q. We would 
however like to remind that the derivation could follow 
other lines of reasoning based on the well known re- 
sults about the evaluation of the entropy of entanglement 

Si,j n of a block of N sites {h,h, ■ ■ ■ ,In}- In fact, 

in Refs. [15|,[l3], it is shown how to express 5^^,..,^ in 
terms of the eigenvalues v\ of a matrix Sn — (TatT^) 1 / 2 . 
There Tjy is the sub-block of order N, relative to the 
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N sites belonging to the block, of a matrix (called T 
in Ref. [15|) equal to ^ T 4> and that is nothing but the 
transpose of the T used in this paper. The case of the 
single site entanglement is then a special case of their 
analysis i.e., N = 1 and, due to the circulant nature 
of T for all i one has that Si = Vi = \Tu\ and Si- — 
-[(1 + 1/0/2] log[(l + 1/^/2] - [(1 - i/<)/2] log[(l - 1/0/2]. 

We now pass to the evaluation of the TDL of Tu for 
L odd. We observe that, since the single particle eigen- 
values are different for j even and odd, we can split 
the sum in T« into an even and an odd part and, since 
in the TDL J^k ~~ * (L/tt) J , we can write these sums as 
integrals such that: 

rp 1 (M-l)dz 



^ 7o v/(^- l) 2 +7 2 cot 2 



r/2 



(/i — l)dx 



7o VO" - I) 2 +7 2 tan 2 



(41) 



The solution of these integral is different inside and out- 
side the region I7I > \\i — 1|. One can write 



2 11 — 1 
T- = — ' 



7^2-^-1)2 



:hl 



V 7 2 - - l) 2 + | 7 | 

Im-i| 



(42) 



for I7I > \n — 1| and 
2 /i- 1 



7T y/{fx— l) 2 — 7' 



• V(/^-l) 2 -7 2 
arcsm — : : (43) 

Ia* — 1| 



for I7I < \n — 1|. T u can be used to study the transition 
H — > 1, while T u will be used to study the transition 
7^0. As far as \x — > 1 is concerned, we first observe 

that T' u (2/7r)ln|/x-l|(/x-l)/|7| -> 0, thus n -► 1/2 

(half filling) and 5^ — > 1 is maximal On the contrary, 
when 7 — > we have that — > sign(/x— 1), consequently 
Si — ► and n tends to zero or one, so that the ground 
state is factorized. 

The next step is the study of the derivatives of 
Tu and Si. One always has: d x Si — (l/2)[log((l — 
|T«|)/2) - log((l - |T«|)/2)]0 B |T«|. We first analyze the 
case x = fi (7 / 0). When /1 — > 1 we have that 
d^rr = -2d^n = d^/L » -ln| M - 1|/| 7 |. The di- 
vergence of the density, that correctly signals the un- 
dergoing phase transition, is however not transferred to 
the first derivative of the single site entanglement; the 
latter has a maximum in fj, = 1 and in fact we have 

dpSt (x log ^gj^l^l w ft - 1) log 2 |M - 1| = 0. The 
divergence is shifted to the derivative of second order, 
i.e., d 2 <W~-log 2 |M-l|- 

We now analyze the case x — 7 (/1 / 1). At the 
transition 7 — » the derivative d 7 T(- — > — sign(7)2[7r(/i — 
1)] _1 is finite but, since when 7 — > 0, T/j' w 1 — 2-f/(ir([i — 
1)) the derivative of the single site entanglement d^Si 
diverges as -sign( 7 ) log (\j\/n\fi - - 1|. 



In summary we can make the following considerations. 
Due to its direct link to the density n, the behavior of 
Si in the TDL reflects the properties of the ground state: 
for 7 = (fi / 1) the state is factorized, while for /i = 1 
(7 =/= 0) the single site quantum correlations are maxi- 
mal. At the transitions the critical behavior is described 
by the derivatives of St. For jj, — * 1 the divergence of 
d 2 Si/dfi 2 is again directly linked to the divergence of 
df/n. In the case 7 — > 0, the derivative 9 7 T/ 4 ' is finite and 
hence it is the functional form of the chosen measure of 
entanglement Si that is responsible for the divergence. 
Indeed d 7 Si oc — logn + log(l — n) and in this transition 
ri — * 0, 1 depending on sign(/i — 1). 



VI. CONCLUSIONS 

In this paper we have given an extensive discussion 
of the relation between ground state fidelity and quan- 
tum phase transitions in quadratic Fermi systems. The 
presented material covers several aspects. First (i) we 
have provided a detailed description of the ground state 
calculation and a simple formula for the fidelity between 
ground states corresponding to different parameters. The 
latter expression is based on the orthogonal part T of 
the polar decomposition of the coupling constant matrix 
Z. The fidelity behavior has then been characterized 
through a combination of its second derivatives, namely 
the minimum eigenvalue h of the Hessian matrix. Subse- 
quently (ii) we have introduced a class of models which 
encompass the XY fermionic Hamiltonian and are an- 
alytically solvable in the cyclic case, thereby offering a 
useful mean to exemplify the above concepts. In particu- 
lar, we have focused on the fully connected system, whose 
long range nature has been shown to give rise to a pe- 
culiar gapful quantum phase transition, where h exhibits 
critical finite-size scaling properties. Finally (iii) for the 
latter model in the cyclic case we have also provided an 
analysis of the single site entanglement, which can be 
extracted from the same matrix T used to calculate the 
fidelity. 

The fidelity approach to QPTs, besides being based on 
an intuitive understanding of the ground state dramatic 
change in critical regions, seems to provide further con- 
ceptual insight into these phenomena, making explicit the 
connections among the ground state, its energy, and the 
single particle spectrum. In particular, the possibility of 
expressing the ground state of these systems only as a 
function of the matrix T completely clarifies its relation 
with the single particle energies. The latter are seen to 
give rise to divergences in fidelity derivatives when some 
of them either tend to zero or to infinite in the thermo- 
dynamic limit. The analytical solvable models presented 
here allow for a study of these relations in a thorough 
way. 

We deem it appropriate to add some remarks on the 
possible developments of this approach. In all the models 
considered so far, a one to one correspondence between 
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QPTs and fidelity drops has been observed. It would 
however be interesting to study the fidelity behavior in 
the presence of more subtle transitions, as in the case of 
topological QPTs, or even in crossover regions, where no 
discontinuous transition is present. In fact, it would be 
desirable to answer the question, whether fidelity drops 
give a reliable signature of QPTs or if they can be orig- 
inated also by to other effects. In other words, could 
fidelity drops (or, more precisely, divergences of fidelity 
derivatives) provide a tool to define QPTs? In spite of 
the possible computational difficulty of the fidelity analy- 
sis in general systems, the added value due to the concep- 
tual simplicity of the ideas underlying our method makes 
the latter, we believe, an interesting question. 



APPENDIX A 

In this appendix we prove the relationship between 
det T and the ground state number parity. 

We start by proving that det T = 1 implies even par- 
ity. If detT = 1 and -1 g SpT then T + 1 is invert- 
ible, G exists, and the even parity follows immediately. 
If instead detT = 1 but — 1 e SpT, then this nega- 
tive eigenvalue must appear an even number of times 
in the spectrum. We note that, if detT = 1, one can 
write T = e , where K is an antisymmetric real ma- 
trix. Then, a real orthogonal transformation R exists 
such that R T KR, and hence R T TR, is block diagonal. 
For L even, one has only 2x2 blocks, while for L odd 
an additional lxl block is present, corresponding to 
the only unpaired eigenvalue, which is for K and 1 for 
T. Being present an even number of —1 eigenvalues in 
SpT, one can group them in pairs in blocks of the form 
— 1 2 . The idea is now to find a parity preserving canon- 
ical transformation of the c$ 's which flips the sign of the 
— 1 eigenvalues of T. In this way one gets a matrix T' s.t. 
detT = 1 and — 1 ^ SpT', being thus possible to define 
G' = (T' - 1)/(T' + 1). The ground state will then have 
the form (J7J with the transformed operators c\ and the 
corresponding vacuum |0'), so that in terms of the new 



parity operator P^i = (—1) , where N' = c 'j ^j: 

the ground state is even. The parity preserving nature of 
the canonical transformation will ensure that Pn 1 = Pn 
and the proposition be proved. 

The transformation properties of the matrix Z for 
a general real canonical transformation have been de- 
scribed in the first part of Sec. [H] We then consider 
the canonical transformation c'j = Rj k Ck, i.e., u = R T , 
v = 0. This transformation is evidently number pre- 
serving, N' = Y%=i c 'j c j — Z)£=i c j°j = N - Further- 
more one has Z' — R T ZR = R T VZZ T RR T TR so that 
T 1 = R T TR is block diagonal. 

Without loss of generality, we now assume that the first 
block of T' is given by —la an d that no other —1 appears 
in the spectrum. We then use a second trivial canonical 



transformation defined by v! + v' — 1 and u' — v' = 
-l 2 ©li_2< In this way one gets Z" = Z'(u'-v') andT" 
is equal to T' apart from the first two diagonal elements, 
whose sign is now changed. Explicitly, u' — ©2 © 1l-2, 
v' = l2ffiCz,_2) an d it is immediate to check that v! and v' 
satisfy the canonical conditions © ■ It is easy to see that 
this transformation, although not number conserving, is 
parity preserving. Indeed one finds N" = N' + 2(1 — n[ — 
n' 2 ), where n'j = c'jCj, so that P^" = Pn> = P^. 

This proves that if det T = 1 the ground state is even. 
Conversely, if det T = — 1 , one can find a trivial parity 
flipping canonical transformation which yields T' such 
that detT' = 1. For example, one can choose ^l| c[ = c\ 
and c'j = Cj for j = 2, . . . ,L, i.e., u — © 1l-i an d 
v = 1 © Ol-i- Then u + v — 1 and u — v = -1 © 
li_i, so that V = T(-l © l L -i) and hence detT' = 
-detT = 1. In addition, N' = N + 1 - 2n x , so that 
the parity is flipped, Pn 1 = —Pn- The ground state will 
then be of the form @ with the transformed operators 
and vacuum |0') = cj|0), being hence of even parity in 
the Cj representation and thus odd in the original one. 
This completes our proof. 



APPENDIX B 



We briefly sketch the main steps of the calculation of 
T'(0) = <9 7 T| 7=0 used in Eq. of Subsec. ITVBl In 
the following, primed quantities imply derivation with 
respect to 7. 

By explicitly deriving T = $ T A _1 $Z with respect to 
7 and using the fact that A'(0) = (D one finds T'(0) = 
[$' T (0)$(0), A^]A - A~ l B', where B' = B/7. 

To avoid the explicit calculation of $'(0) one can use 
the relation [A,B'] = [$' T (0)$(0), A 2 ], obtained from 
(ZZ T )'(0) with the aid of A'(0) = 0. By noting that 
A 2 = [L + 2(>-l)]A-(L + M-l)(M-l)l [seeEq. Jg)] 
and that [$' T (0)$(0), A^]A = -A^ 1 [$' T (0)$(0), A] 
one finds 



T'(0) 



-A' 



[A,B'J 



L + 2(n-l) 



B' 



(Bl) 



which, after algebraic manipulations, yields Tj fc (0) = 
[(fc-j)/(i/2+ M -l)-sign(fc-j)]/(M-l). Then, Eq. ® 
is obtained by calculating only the diagonal elements of 
[T'(0)] 2 , namely 



1 



2(A* - i)iC? -L-l) 



(L/2 + ix - 
l T(T + l)(T + 6^- 2) 
6 (L/2 + fi- l) 2 



I) 2 
-L 



(B2) 



and finally summing to get the trace. 



14 



[1] A. Osterloh, L. Amico, G. Falci, and R. Fazio, Nature 
416, 608 (2002); J. Vidal, G. Palacios, and R. Mossery, 
Phys. Rev. A 69, 022107 (2004); J. Vidal, R. Mossery, 
and J. Dukelsky Phys. Rev. A 69, 054101 (2004); 
M.M. Wol f, G. O rtiz, F. Verstraete, and J.I. Cirac, 
|cond-mat /0512180 

[21 Shi-Liang Zhu Phys. Rev. Lett 96, 077206 (2006); 
A. Hamma, quant-ph/0602091 J.K. Pachos and 
A. CM. Carollo lquant-ph/0602i54| 

[3] S. Ryu and Y. Hatsugai, cond-mat/0601237 

[4] H.T. Quan et a l, Phys. Rev. Lett. 96, 140604 (2006); 
D. Rossini et al, quant-ph/0605051 

[5] E. Farhi et al, Science 292 472 (2001). 

[6] A. Hamma and D.A. Lidar, quant-ph/0607145 

[7] P. Zanardi and N. Paunkovic, |quant-ph/0512249| Phys. 
Rev. E, to be published. 

[8] P. Zanardi, M. Cozzini, and P. Giorda, 
quant-ph/0606130 

[9] E. Lieb, T. Schulz, and D. Mattis, Ann. of Phys. 16, 407 
(1961). 

[10] D. Rossini et al, cond-mat/0605051 

[11] M.-C. Chung and 1. Peschel, Phys. Rev. B 64, 064412 

(2001); I. Peschel, J. Phys. A 36, L205 (2003). 
[12] R. Bathia, Matrix Analysis, Springer- Verlag, New York, 

Inc., 1997. 

[13] A. Perelemov, Generalized coherent states and their ap- 
plications, Springer- Verlag , New York, Inc., 1986. 

[14] E. Barouch, B.M. McCoy, and M. Dresden, Phys. Rev. A 
2, 1075 (1970); E. Barouch and B.M. McCoy, Phys. Rev. 
A 3, 786 (1971). 

[15] J. P. Keating and F. Mezzadri, Commun. Math. Phys. 
252, 543 (2004). 

[16] H.J. Lipkin, N. Meshkov, and A.J. Click, Nucl. Phys. 62, 
188 (1965). 

[17] G. Vidal, J.I. Latorre, E. Rico, and A. Kitaev, Phys. Rev. 
Lett. 90, (2003) 227902. 



[18] A. Anfossi, P. Giorda, A. Montorsi, and F. Traversa, 
Phys. Rev. Lett. 95, 056402 (2005); A. Anfossi, C.D.E. 
Boschi, A. Montorsi, and F. Ortolani, Phys. Rev. B 73, 
085113 (2006); S. Gu et al., Phys. Rev. Lett. 93, 086402 
(2004); D. Larsson and H. Johannesson, Phys. Rev. A 
73, 042320 (2006); R. Somma, G. Ortiz, H. Barnum, E. 
Knill, and L. Viola, Phys. Rev. A 70, 042311 (2004). 

[19] P. Zanardi, X. Wang, J. Phys. A 35, 7947 (2002); P. 
Zanardi, Phys. Rev. A 65, 042101 (2002). 

[20] Recalling that the determinant of a matrix only changes 
its sign under the exchange of matrix rows or columns, 
one can rearrange the four L x L blocks which build C — 
CI21, so to prove the implication. 

[21] This is due to the anticommutation rules of fermionic op- 
erators: for L even c|c^ , 2+1 must have the same weight 
as <^ l / 2 +i c l+i = c' L , 2+1 c\, because of periodicity, but 
this in turn is nothing but ~ c |cj^ 2+1 , and the two terms 
cancel in the Hamiltonian. Then, the corresponding en- 
tries in the antisymmetric matrix B have to vanish. 

[22] Recall that /i = 2A in terms of the Hamiltonian 1161 . 

[23] At (fi = 0,7 = 1) the matrix ZZ T can be diago- 
nalized analytically for any L. The eigenvalues, apart 
from the lowest one which is always zero, obey the 
formula 1 + tan 2 [jV/(2L - 1)], j = 1,...,L - 1. In 
the standard basis, the components of the normalized 
eigenvector corresponding to the j-th non zero eigen- 
value are 2sin[j7r(2fc - 1)/(2L - 1)] cos(kir)/y/2L - 1 for 
k = 1, . . . , L — 1 and zero for k — L. The only non-zero 
component of the zero eigenvalue eigenvector is instead 
the L-th one. 

[24] Similar first order QPTs happen also in the usual XY- 
model, as for example at /i = —2 (i.e., A = — 1) for L odd 
in the c-cyclic case. 



